q-deformed logistic map with delay feedback 
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The delay logistic map with two types of q-deformations: Tsallis and Quantum-group type are 
studied. The stability of the map and its bifurcation scheme is analyzed as a function of the 
deformation and delay feedback parameters. Chaos is suppressed in a certain region of deformation 
and feedback parameter space. The steady state obtained by delay feedback is maintained in one 
type of deformation while chaotic behavior is recovered in another type with increasing delay. 



I. INTRODUCTION 

Theory of quantum integrable systems [J, [2[ has ini- 
tiated a new type of symmetry and associated with it 
mathematical objects called quantum groups. These are 
related to the usual Lie groups as quantum mechanics is 
related to its classical limit. Physically a group quanti- 
zation can lead to a kind of deformation of the group 
manifold, related to a physical (classical or quantum) 
system, q-deformation of many classical Lie groups has 
stimulated much activity in the pursuit of understand- 
ing its physical meaning, due to the emergence of quan- 
tum group like features in many physical systems. It has 
been realized that q-deformation effectively takes into 
account the interactions in physical systems 3-6]. The 
q-deformation is non-trivial in the sense that the emerg- 
ing deformed algebra is no longer linear. It would be 
constructive to study q-deformation in the context of 
dynamical systems. 

One of the popular model of discrete nonlinear dynam- 
ical systems is logistic map [7J . The study of dynamical 
system with delay is important when the information is 
feedback with a measurable delay, e.g., due to the spa- 
tial extensiveness of the system and the finite velocity 
of propagation of information, or when the characteristic 
timescale of the system is smaller than the delay time Q . 
The delayed equations have been used for modeling pur- 
poses in optics 0, [l(| , chemistry [llj , and biological sys- 
tems [lJjLUl. The delay feedback also plays an important 
role in controlling chaos [1J, [l5[ . The interplay of delay 
and nonlinearity plays a central role in self-organization 
and complex phenomena of dynamical systems and chaos 
is suppressed or controlled by stabilizing unstable peri- 
odic orbits with delay feedback [16—20]. In another work 
[2l| , the normal logistic and exponential maps were used 
to study the transition from chaotic to regular dynamics 
induced by stochastic driving. 

A one-dimensional logistic map is a non-linear differ- 
ence equation 



kn+l 



O^X n \l X n ) : 



(1) 



where a is a constant, and is taken to be positive in the 
rest of the paper. Also, x n denotes the value of x after n 
iterations. Eq. (JT|) arises, for e.g., in the case of modeling 



of population growth 



dN(t) 
dt 



r(t)N(t), where N(t) is the 



population at a time t and r(t) is the difference between 
birth and death rates per head of the population. 

As q-deformation essentially involves modification of 
a function such that in the limit of q — > 1 the usual func- 
tion is obtained, there is no unique q-deformation for a 
function. On the other hand, delay transforms the dy- 
namical state of the system. It is therefore natural to 
use q-deformations suitably in the study of non-linear 
systems with delay feedback. Here we discuss two forms 
of q-deformations of the logistic map, studied in the lit- 
erature, but with the additional proviso that the map has 
a memory inbuilt into it, in the form of a feedback mech- 
anism. This has the advantage of studying the system 
from a more realistic perspective as well as the possibil- 
ity of having a chaos suppression mechanism inbuilt into 
the system. 

The paper is setup as follows. In Section II, we discuss 
the basic deformations, that will be studied, along with 
delay. The logistic equation deformed by the two pre- 
scriptions and undergoing a delayed feedback are studied 
in Sections III (A) and (B), respectively. Along with 
a numerical study of the interplay between the various 
parameters, a stability study of the two systems is also 
made. Section IV concludes the paper. 



II. g-DEFORMATION WITH DELAY 



A deformation of the logistic map (E 
the non-extensive statistics of Tsallis [2! 
[23J J . in which the map 

x n+ i = a[x„](l - [x„]), 



. Q}, based on 
was proposed 



was considered. Here [x] = 
2 for x in the interval [0, 1 



(2) 

TT7i — %i — r> and — oo < q < 
l+(l— g)(l— x) ' » 

An important difference be- 
tween |T]) and j2| is that the deformed map (0 is concave 
in parts of x-space while the map without deformation 
(U|) is always convex. Further, the use of @ showed the 
rare phenomena of the co-existence of attractors, i.e., the 
co-existence of normal and chaotic behavior. 

The logistic map with delay (r) feedback is given by: 



I-ra+l 



= f(x n ) = aor„(l - x n )(l - /3) + f3x n - T (3) 




FIG. 1: The parameter space (e, /3) showing chaotic (black 
color), periodic (brown color for period > 2 and blue for pe- 
riod 2), period 1 (non-zero steady state in blue color) and 
period 1 (zero steady state in white). Here nonlinearity pa- 
rameter a — 4, delay parameter r = 1. As seen from the 
figure, with increase in coupling (/?), for fixed q-deformation 
e, the system undergoes a transition from the chaotic to the 
periodic regime, thereby highlighting the role played by the 
feedback on chaos suppression. 



where, /? is the feedback amplitude and r is the delay 
time. The analogous q-deformed map with delay is: 



X n +1 = a[x n ](l ~ [Xn])0- ~ P) + P[x n - 



(4) 



Another Quantum-group (Qu-group) type of q- 
deformation, of the logistic map, was proposed [24J as: 

(5) 



[x n +i] = a[x„](l - [x n ]), 



where 



(6) 



Here q is real and x is in the interval [0,1]. This q- 
deformed logistic map is different from Eq. ([2]). It is not 
possible to transform the q-deformed map, introduced in 
Eq. ([5]), to that in Eq. ([2J. In particular, it is not possible 
to relate the q parameter in Eq. ([2]) to the q parameter 
in Eq. ([5]). Further, in the proposed q-deformed map 
Eq. ©, the left hand side is also q-deformed, in contrast 
to Eq. ([2]). Thus in the space of q-deformed variables 
the q-deformed logistic map Eq. ([5]) , is the usual map. In 
mapping to ordinary space, all the corresponding phys- 
ical features emerge. This is not possible in the map 




FIG. 2: a) Bifurcation diagram depicting x with respect to fi 
for a = 4, t = 1. As seen from the figure, with increase in 
coupling (f3), for fixed q-deformation e — 0.5, the system un- 
dergoes a transition from the chaotic to the periodic regime, 
ultimately to the steady state, confirming the scenario de- 
picted in Fig. © b), where the Lyapunov exponent A is 
plotted with respect to j3 for a = 4, r = 1. As seen from the 
figure, with increase in coupling (/?), for fixed q-deformation 
e, the system undergoes a transition from the chaotic (A > 0) 
to the periodic (A < 0) regime, ultimately to the steady state, 
c) Time series showing chaotic, period-2 and steady states for 
P = 0.0, 0.5, 0.7, and 0.8 respectively. 



Eq. @. In the limit q — > 1, it is seen that [x] — > x and 
we obtain the usual logistic map ([1]). 

The corresponding deformed map with delay would be: 



[X n +l] = «[*n](l - [»n])(l ~ P) + P[x v 



(7) 



What do we expect from such a study? q-deformations 
simulate correlations in the system while the delayed 
feedback brings in memory. An interplay of these two 
effects should help in understanding the mechanism of 
suppression of chaos in systems that have inbuilt corre- 
lations. 



III. ANALYSIS OF THE g-DEFORMED 
LOGISTIC MAP WITH DELAY 

Here we take up the logistic map with delayed feedback 
and q-deformed according to both the prescriptions, i.e. 
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FIG. 3: Bifurcation diagram depicting x with respect to (3 for 
a — 4 q — 0.5, r = 1, 2 and 3. System behavior is period 
two cycle, fixed point and chaotic for higher delay feedback 
strength for r = 1, 2 and 3 respectively. In the inset, is plotted 
the critical value of the feedback parameter /3 C , for which the 
system goes into the steady state x* — 0. 



according to Eqs. ([2]) (Tsallis type) and (0 (Qu-group 
type), respectively. 



A. Tsallis type of deformation 

1. Stability Analysis: Analytical Results 

We make an analytical study of the effect of memory 
in Eq. (J4|) . This involves expanding the original equation 
to a set of r + 1 equations, r being the delay [8(. In order 
to capture the essence of our q-deformed delayed logistic 
map, we take up the case of r = 1. For the one-cycle 
stability, using Eq. (|4]): 



X n+X = a(l-0)(2-q) 



&n\*- •^n) 



[l + (l-g)(l -*„)]* 



-P 



Vn 



l + (l- q )(l-y n )- 



y n+1 = x n . (8) 

The corresponding Jacobian matrix takes the form: 



J(x,y) 



2 — q—(3 — q)x n 



m-q) 



From the characteristic equation of Eq. (fTU|) . its eigen- 
values are: 

\ x 2 = a{ ) ~ P] ± l Ja 2 (l-f3) 2 +4p(2-q). 
2(2 -q) 2(2 -q) V y H > Hy H ' 

(11) 
From the above eigenvalues, the condition for the stabil- 
ity of the fixed point is obtained as: 



= a{i p)(z q) [1+(1 _ g)(1 _ a . n)] 3 [i+(i_ g )(i_3, n )] 



which for the trivial fixed point: x — y = gives: 

/ a(l-/3) g \ 



(9) 



(3> 



a - (2 - g) 
a-1 



(12) 



These results are borne out by the numerical results 
shown below, where along with one-cycle stability, multi- 
cycle stability is also analyzed. 



2. Numerical results 

We fixed the nonlinearity parameter a = 4 and delay 
time r = 1 to study the effect of q-deformation and delay 
feedback in (q, (3) parameter space as shown in Fig.[TJ For 
lower feedback strength j3, system is in chaotic (in black) 
or in higher period (in brown) state. With sufficiently 
large delay feedback, system goes to period two cycle 
(in blue) and then to steady state x* ^ (in red) and 
then to x* — (in white). The bifurcation diagram and 
Lyapunov exponent as a function of feedback strength 
j3 is plotted in Figs. HJb) and (c), respectively, for fixed 
q-deformation q — 0.5 which confirms chaos suppression 
via reverse period-doubling bifurcation. Time series in 
different dynamical states, i.e. chaotic, period two cycle, 
and steady states are shown in Fig.[^c). When the delay 
time r is further increased, the transition from chaos to 
steady state is seen to occur at lower values of feedback 
strength j3 and the critical value of feedback strength d c , 
when the system approaches the steady state x* — 0, is 
0.76 for t — \ and approaches to 0.63 with increasing r 
(see Fig. H. 

Thus, in the Tsallis- type of q-deformed logistic map 
with delay, with increase in feedback, stability is seen 
for all delays; in contrast to the situation in the corre- 
sponding map without deformation, where stability was 
observed for only odd delays [8j. 



B. Quantum— group type of deformation 

1. Stability Analysis: Analytical Results 



.7(0,0) 



(2-q) (2-q) 

1 



We make an analytical study of the effect of memory 

in Eq. ([7]). As before, in order to capture the essence of 

qq~) our q-deformed delayed logistic map, we take up the case 

of r = 1. For the one-cycle stability, using Eq. (J7J|, we 
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FIG. 4: The parameter space (q, f3) showing chaotic (black 
color) and periodic (brown for period > 2 and gray for pe- 
riod 2). Here nonlinearity parameter a = 4, delay parameter 
T — 1. As seen from the figure, with increase in coupling 
(/3), for any given q-deformation e, the system undergoes a 
transition from the chaotic to the periodic-2 cycle regime, 
thereby highlighting the role played by the feedback on chaos 
suppression. 



get the following two coupled equations: 
a(l - (3) 



X n +i — 1 



(l-Q) 



(l-X n )(X n -q)-f3(l-Y n ) 

(13) 



where, for convenience in calculations, we have made the 
change in variable : q Xn — X n . The corresponding Jaco- 
bian matrix takes the form: 

j {x , Y) =(f^[2x„-( q+ i)]e\ (14) 

which for the trivial fixed point: x = y — 0, correspond- 
ing in the new variables to: X = Y = I, gives: 



.7(0,0) 



a(l-/3) p 
1 



(15) 



From the characteristic equation of Eq. (fT5|) . its eigen- 
values are: 

Ax,a = " (1 ~ P) ± ~ V« 2 (l ~ P) 2 + 4,3. (16) 
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FIG. 5: a) Bifurcation diagram depicting x with respect to 
p for a — 4, r = 1. As seen from the figure, with increase 
in coupling (/3), for fixed q-deformation q = 0.5, the system 
undergoes a transition from the chaotic to the periodic-2 cy- 
cle regime confirming the scenario depicted in Fig. Q. b) 
Lyapunov exponent A with respect to /3 for a — 4, t = 1. As 
seen from the figure, with increase in coupling (/}), for fixed 
q-deformation q, the system undergoes a transition from the 
chaotic (A > 0) to the periodic (A < 0) regime, ultimately to 
the period-2 cycle state, c) Time series showing chaotic and 
period-2 states for /3 — 0.0 and 0.4 respectively. 



From the above eigenvalues, the condition for the stabil- 
ity of the fixed point is obtained as: 



/3<1. 



(17) 



Since these calculations are made in the new variable X n , 
which is related to the original variable x n by x n = " " ; 
the result of Eq. (|T71) when interpreted in the original 
variable implies that x* — will never stabilize under 
positive feedback. This is borne out by the numerical 
results shown below, where the one-cycle as well as multi- 
cycle stability is analyzed. 



2. Numerical results 

Again, we fixed the nonlinearity parameter a = 4 and 
delay time r = 1 to study the effect of q-deformation 
and delay feedback in (q, j3) parameter space, as shown 
in Fig. 2J For lower feedback strength /3, system is in 




back strength and with further increase of delay time 
system goes back to chaotic state as shown in Fig. [6) For 
t = 2, within period-1 cycle, some values of feedback 
strength (3 take the system to period-3 cycle as shown 
in Fig. [6th). It shows that chaos cannot be suppressed 
in Qu-group type of deformation for delay feedback with 
r > 2. We have studied the asymmetric Qu-group type 
of deformation case (x n +i = f([x n ])) also and find simi- 
lar results (not shown here). Thus here the non-linearity 
is predominantly due to the form of the deformation cho- 
sen. In contrast to Tsallis type of deformation, the Qu- 
group type of deformation, suppress chaos to two-cycle 
as compared to one-cycle. 



IV. CONCLUSIONS 



FIG. 6: Bifurcation diagram depicting x with respect to j3 for 
a — 4 and q = 0.5, for increasing delay time: r = 1, 2 and 3, 
respectively, a) Period-2 cycle for r = 1, b) fixed point for 
t = 2, and c) chaotic behavior for r = 3 are observed. 



chaotic (in black) or in higher period (in brown) state. 
With sufficiently large delay feedback, system goes to 
period two cycle (in gray). The bifurcation diagram and 
Lyapunov exponent as a function of feedback strength (3 
is plotted in Figs. [5ja) and (b), respectively for fixed q- 
deformation q = 0.5 which confirms chaos suppression via 
reverse period-doubling bifurcation to period two-cycle. 
Time series in different dynamical states, i.e., chaotic and 
period two cycle are shown in Fig. EJc). For delay time 
t = 2, system goes to steady state with sufficient feed- 



In this work, we have studied the logistic map from the 
perspective of q-deformation and delay feedback. This 
enables us to study the interplay between the competitive 
features, viz. complexity, brought by the q-deformation 
and order, by delay feedback. Chaos is suppressed with 
feedback for both kinds of deformations, i.e., Tsallis and 
Qu-group type. However, for the Tsallis type of defor- 
mation, a steady state is achieved, an observation which 
validates its use in statistical mechanics of complex sys- 
tems, where one would expect a system to eventually go 
to a steady state, while for the Qu-group type of defor- 
mation the period-2 cycle gets stabilized. With increas- 
ing delay time, the transition to steady state is obtained 
at lower values of feedback strength in Tsallis type of q- 
deformation while chaotic state is recovered in Qu-group 
type of q-deformation with increasing delay time. 



[1] V. Drinfeld Dokl. Akad. Nauk. 283 1060 (1985). 
[2] M. Jimbo Lett. Math. Phys. 10 63 (1985). [15 

[3] G. J. Estere, C. Tejel and B. E. Villarroya J. Chem. Phys. [16] M 
96 5614 (1992). [17] M 

[4] Z. Chang and H. Yan Phys. Lett. A 154 254 (1991). [18] 

[5] Z. Chang, H. Y. Guo and H. Yan Phys. Lett. A 156 192 [19] 
(1991). [20] 

[6] Parthasarathy R 1993 LMSc-93/23, Preprint. 
[7] Jose J V and Saletan E J 2002 Classical Dynamics: A [21] 

Contemporary Approach (Cambridge University Press) 
[8] T. Buchner et al. , Phys. Rev. E 63, 016210 (2000). [22] 

[9] K. Ikeda, Opt. Commun, 30, 257 (1979). 
[10] G. Giacomelli, R. Meucci, A. Politi, and F. T. Arecchi, 

Phys. Rev. Lett. 73, 1099 (1994). 
[11] I. R. Epstein, J. Chem. Phys., 92, 1702 (1990). [23] 

[12] M. C. Mackey and L. Glass, Science, 197, 287 (1977). 
[13] S. Cavalcanti and E. Belardinelli, IEEE Trans. Biomed. [24] 

Eng., 43, 982 (1996). 
[14] E. Ott, C. Grebogi, and Y. A. Yorke, Phys. Rev. Lett. 



64, 1196 (1990). 

K. Pyragas, Phys. Lett. A 170, 421 (1992). 

de Sousa Vieira et al. , Phys. Rev. E 54, 1200 (1996). 
M. D. Shrimali, et al. , Phys. Lett. A 374, 2636 (2010). 
E. Fick et al. , Phys. Rev. A 44, 2469 (1991). 
C. Masoller et al. , Phil. Trans. R. Soc. 369, 425, (2011). 
Complex Time-Delay Systems. F. M. Atay (ed.), 
Springer- Verlag Berlin, 2010. 

Prasad A and Ramaswamy R 1999 

eprint\arXiv:chao-dyn/99 1 1002 

C Tsallis 1988, J. Stat. Phys. 52, 479; M Gell-Mann 
and C Tsallis (Eds.), Nonextensive entropy - Interdisci- 
plinary applications (Oxford University Press, New York, 
2004) 

R. Jaganathan and S. Sinha, Phys. Lett. A 338 277 
(2005). 

S. Banerjee and R. Parthasarathy J. Phys. A: Math. 
Theor. 44 045104 (2011). 



